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We analyze theoretically complete wetting of a substrate supporting an array of parallel, vertical 
plates which can tilt elastically. The adsorbed liquid tilts the plates, inducing clustering, and thus 
modifies the substrate geometry. In turn, this change in geometry alters the wetting properties of 
the substrate and, consequently, the adsorption of liquid. This geometry-wetting feedback loop leads 
to stepped adsorption isotherms with each step corresponding to an abrupt change in the substrate 
geometry. We discuss how this can be used for constructing substrates with tunable wetting and 
adsorption properties. 



I. INTRODUCTION 

The scientific interest in wetting properties of liquids 
at solid surfaces dates back to the 18th century. Sev- 
eral reviews account for the progress in understanding 
such wetting phenomena [H-Q . Recent important devel- 
opments include the wetting of rough surfaces and su- 
perhydrophobicity d, Q and the attempts to build "lab- 
on-a-chip" devices, in which liquids are controlled at the 
micro- and nanoscale 0, [1| ■ 

A large part of this recent development is concerned 
with designing surfaces in order to obtain specific, tun- 
able wetting properties, leading to "smart surfaces" . In 
this spirit non-rigid substrates literally provide additional 
flexibility in tailoring wetting properties. Recently the 
usefulness of the interplay of surface elasticity and wet- 
ting has been demonstrated in various accounts @ . They 
include capillary-driven self-assembly [l0l - [l6j . enhanced 
condensation properties , control of surface color [l^ , 
measurement of pressure in nano-channcls [l9| . and a 
biomimet ic p roposal to achieve superhydrophobic sur- 
faces [U . This last suggestion has been refuted [l^l 
but calculations indicate that a related effect might work 
(23| . Even at the macroscopic level the interplay of wet- 
ting and elasticity gives rise to rich phenomena [23 - [26| . 

Here we study an elastic surface structure which 
changes its geometry upon adsorption. In turn, this 
change in the topography of the substrate affects the ad- 
sorption properties. Our goal is to give a bird's eye view 
of the influence of elasticity on wetting and adsorption 
phenomena rather than to calculate specific, detailed re- 
sults. As such, we choose one of the simplest possible ex- 
amples of the influence of elasticity: a surface decorated 
with rigid, vertical, planar plates which can pivot elasti- 
cally upon their base. This choice simplifies the analysis, 
while still retaining the important elements of the phe- 
nomenology. The adsorption isotherms exhibit jumps at 
undersaturations at which the surface changes its geom- 
etry. The different kinds of possible isotherms depend on 



elasticity, geometry, and surface tension. We discuss how 
the use of this kind of surface together with the applica- 
tion of an electric field can lead to "smart surfaces" with 
"on-the-fly" tunable wetting properties. 

In the next section we introduce the phenomenology 
of wetting, filling, and capillary condensation and qual- 
itatively describe the wetting of flexible substrates. In 
Sect, mil we present the results for our simple macro- 
scopic model. We close with discussions of the results 
(Sect. lIVp and of their implementation in order to ob- 
tain surfaces with tunable properties (Sect. |V| . 



II. PHENOMENOLOGY 

Before we consider the wetting and adsorption prop- 
erties of elastic substrates we first briefly review some 
of the most important results for wetting of planar and 
patterned substrates (such as wedge flUing and capillary 
condensation), because we rely heavily on these for our 
calculations. 

In all situations we consider a solid substrate in ther- 
mal equilibrium with a bulk gas phase at a certain tem- 
perature T and pressure P (or chemical potential /i), thus 
considering the grand canonical ensemble. Away from 
gas-liquid coexistence in the bulk a thin layer of liquid 
adsorbs on the substrate. The adsorption isotherm is 
obtained by monitoring the excess adsorption F upon in- 
creasing the pressure towards liquid-vapor coexistence in 
the bulk at fixed temperature. The excess adsorption 
can be measured experimentally with hi gh accuracy by 
using, e.g., microbalance techniques |27H31| . This is a 
classical and one of the easiest methods to characterize 
the wetting properties of a substrate. In the present con- 
text we consider the number densities pi and py of the 
liquid and vapor phase, respectively, to be spatially ho- 
mogeneous so that the excess adsorption equals pi — py 
times the volume occupied by the liquid phase. Factoring 
out Ap ^ pi ^ Py and the linear extension of the system 
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FIG. 1: Characteristic adsorption isotherms of a planar sub- 
strate as the cliemical potential ^ approaches bulk gas-liquid 
coexistence /io(r) at constant temperature. F is the excess 
amount of liquid adsorbed for partial wetting (a) and com- 
plete wetting (b). In (b) the type of divergence of F reflects 
the asymptotic decay of the underlying interactions in the 
system [3. 

in the translationally invariant direction, F is the area 
occupied by liquid in the cross section (see Appendix A) . 

For a flat substrate there are two typical behaviors of 
the adsorption isotherms (see Fig. [T]): either the excess 
adsorption F remains finite up to gas-liquid coexistence 
fio{T) or it diverges as coexistence is approached [1]. The 
former behavior is called partial wetting, corresponding 
to a contact angle > 0° at gas-liquid coexistence, and 
the latter is called complete wetting, corresponding to 
9 = 0°. The character of the divergence depends on 
the microscopic details of the molecular interactions. For 
non-retarded van der Waals interactions the adsorption 
diverges as F oc (A/i)""^ [sj, where Afi = fJ,o{T) — fJ, > 
is the deviation of the chemical potential of the bulk gas 
phase from its value fjLo{T) at liquid- vapor coexistence in 
the bulk. 

Due to confinement, fluid will condense in a pore at a 
lower value of the chemical potential, compared to bulk 
coexistence. For the simplest possible pore, formed by 
two planar, parallel, inflnite plates a distance L apart, 
the comparison of the free energies for the empty and the 
filled pore leads to the Kelvin equation, which provides 
the value of the chemical potential at which capillary 
condensation occurs 0] : 

. 2a cos 9 

Am-Am--^^. (1) 

Here a is the surface tension of the liquid- vapor interface. 

Another important element of our later analysis is the 
adsorption at a wedged substrate with a tilt angle a (see 

t.[2Ia)). Macroscopic thermodynamic arguments [32|, 
show that upon raising temperature the filling of a 
wedge precedes wetting of its walls and occurs when the 
contact angle at the planar wall equals the tilt angle of 
the wedge, i.e., 9 = a. Insight into the wetting behavior 
of more complicated substrate geometries can be gained 
by using a simple mesoscopic model of adsorption [sj] . 

With these prerequisites, we are now in a position to 
qualitatively describe the adsorption of liquid at an elas- 



tic substrate. To be specific, we consider a substrate 
decorated with a periodic array of planar plates which 
are vertical if there is no liquid adsorbed (see Fig. 
The underlying substrate is rigid as well as the plates as 
such. However, the plates can tilt elastically around their 
contact line with the underlying substrate. 

We consider also this substrate to be in thermal equi- 
librium with a gas phase at a certain temperature T and 
chemical potential fj. characterizing the grand canonical 
ensemble. If the substrate is stiff (i.e., if elasticity can be 
ignored) the excess adsorption isotherm F typically looks 
like the dotted line in Fig. For large undersaturations 
A/i only a small amount of liquid adsorbs at the walls, 
forming a microscopic wetting film. The precise amount 
of that adsorbed liquid depends on the microscopic de- 
tails of the system and can be calculated with the same 
methods used for a planar substrate At A^ = A/i'^'^ 
the adsorption isotherm exhibits a steep increase at an 
undersaturation corresponding to the onset of capillary 
condensation of liquid in the space between the plates. 
Unlike capillary condensation between parallel walls, this 
increase is actually smooth, albeit very steep, because the 
pore is a capped capillary [35l - l38| . For even smaller un- 
dersaturations the liquid continues to adsorb on the to- 
pographically structured substrate with the type of the 
corresponding divergence being controlled by the range 
of the molecular interactions involved, similar to a pla- 
nar substrate [|. For A^ < An"" and A/i ^ one can 
identify several filling and wetting regimes (35l . l36j . 

We now consider elasticity, i.e., the plates can bend 
due to the capillary forces of the adsorbed liquid. Start- 
ing from large undersaturations, i.e., Afi :S> Afi"", upon 
decreasing A/j, the adsorption isotherms for the elastic 
substrate are de facto the same as the ones for the cor- 
responding rigid structured substrate until capillary con- 
densation occurs. (Concerning flexible slit pores which, 
however, maintain their geometry, see Refs. (39j and [40l|.) 
At the undersaturation Aji"" corresponding to capillary 
condensation the liquid starts to fill the space between 
the plates as before. However, if the elastic constant 
for tilting a plate around its baseline is not too big, the 
ensuing capillary forces can cause the plates to cluster 
together. This can be rationalized by noting that a small 
deviation of a plate from its vertical orientation causes 
the pore to be smaller on one side and bigger on the 
other side; accordingly, these pores form capped wedges 
(see Fig. [IKb)) with tilt angles a smaller or larger than 




FIG. 3: (Color online) Model elastic substrate. The rigid underlying substrate is decorated with rigid, planar, and equally 
spaced plates. The substrate and the plates are considered to be formed by the same material. The plates can pivot elastically 
upon their line-like bottom, where they are anchored at the underlying substrate, represented by a dark gray spot. An 
adsorption-induced collapse of plates triggers a chain of clustering events. The arrows represent the concomitant movement 
of the liquid-vapor interface and of the plates. Subsequently the clustered plates can collapse into bigger clusters. The sketch 
exaggerates the thickness of the plates, which are considered to be thin enough, compared to L, so that the thickness plays no 
role in our calculations but thick enough so that microscopic details, such as interactions between portions of liquid adsorbed 
on different sides of a plate, are not important. 



90°, respectively. In line with the behavior of wedge fill- 
ing (HjUlillH liquid condenses in the smaller pores while 
it evaporates from the larger ones, relative to the upright 
configuration a ~ 90°. The condensation and evapora- 
tion of liquid in pores of different sizes create an imbal- 
ance of forces acting on their common plate due to the 
surface energies of the interfaces involved. This imbal- 
ance triggers a clustering cascade, in which the plates 
form clusters of two plates each (Fig. [3]), assuming that 
the plates are not too rigid. This first plate clustering 
event occurs at the same undersaturation as the afore- 
mentioned capillary condensation: A/i'*'^ = A/i"^^. 

At the end of this process the substrate consists of 
capped wedges, with the plates clustered two-by-two. 
Between each capped wedge there is a filled wedge up- 
side down (Fig. [S]). Upon decreasing the undersatura- 
tion further, adsorption proceeds up to A/x = A^'^^ 
where the filling of these capped wedges causes pairing 
in clusters of four plates by the same instability mecha- 
nism as described above. In line with experimental evi- 
dence p^ - ii^ l we assume that, once formed, clusters can- 
not separate. In this case and at this stage the system 
consists of capped wedges with a smaller tilt angle a, 
which will in turn collapse into clusters of eight plates 
(at A^ = A^(^^), etc. In our notation, upon decreasing 
A/Lt there is a first-order transition between a cluster of 
order i (i.e., containing 2' plates) to a cluster of order 
i + 1 containing 2*+^ plates. The upright configuration 
corresponds to i = 0. 

According to the above description the adsorption 
properties change upon every clustering event, turning a 
smooth adsorption curve into a series of adsorption steps 
which coincide with the abrupt changes of the substrate 



geometry, as sketched in Figdl These collapses are lim- 
ited by the length of the plates or by the elastic energy 
associated with the tilting of the plates, i.e., finally the 
clusters become too rigid. When no further collapses are 
possible the adsorption proceeds by filling of the capped 
wedges (see the smooth increase of F in Fig. [4] below 
the last step) until the liquid completely fills the space 
between the plates. Therefore the value of F at the bot- 
tom end of the last step (corresponding to a filled capped 
wedge configuration) equals the one at the upper end of 
the steeply increasing dotted line corresponding to a filled 
upright configuration. Further adsorption causes the liq- 
uid to submerge the plates and the surface tension of the 
liquid- vapor interface no longer acts on the plates. The 
plates can remain stuck to each other if van der Waals or 
other forces keep them together. Alternatively, elasticity 
can be strong enough to restore the plates to their up- 
right configuration. In the former case further adsorption 
is similar to adsorption on a substrate formed by capped 
wedges, whereas in the latter case further adsorption is 
the same as on the rigid substrate, with the plates in the 
upright configuration; this is the case drawn in Fig. 3] In 
either case, for the slope of this part of the adsorption 
isotherm elasticity does no longer play a role. 

If the substrate is not perfectly periodic but there is 
some randomness in the spacing of the plates the ad- 
sorption will proceed in a less orderly fashion [l^. For 
example, there can be clusters of three plates, or clusters 
of different sizes within the overall configuration. This 
will result in more steps in the adsorption isotherm. 

In the next section we present results for a specific 
model system, substantiating the previous qualitative 
reasoning. 
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FIG. 4: (Color online) Sketch of the excess adsorption 
isotherms F of a rigid, patterned substrate with the plates 
fixed in upright position (dotted red line) and of a flexible 
substrate (continuous black line) with the same geometry if 
there is no adsorption (i.e., A/i ^ Afi""). The latter exhibits 
jumps when the substrate geometry changes, with the order i 
of the clusters indicated. When the liquid completely covers 
the plates of the flexible substrate, the plates can return to 
their original upright position and the adsorption isotherm is 
the same for A/x — >■ 0, or they can remain clustered and the 
isotherm will be different. The plot corresponds to the for- 
mer case. The substrate material exhibits complete wetting 
by the liquid (i.e., 9 = 0). See the main text for a detailed 
discussion. 



III. MODEL CALCULATIONS 

In order to underpin the previous phenonienological 
description, in this section we present calculations for a 
simple model which still captures the most important fea- 
tures but avoids the complexity of more realistic models. 

As stated before our choice of the substrate is idealized 
in that it consists of stiff planar plates which tilt around 
their line of contact with the underlying substrate. A 
more realistic system would consist of elastic plates which 
can adopt also nonplanar shapes. This would lead to 
much more difficult calculations without adding essential 
new features. Moreover, the system we analyze is transla- 
tionally invariant along one spatial direction, thus being 
effectively two-dimensional. This invariance is broken by 
studying systems with arrays of elastic posts. Such an 
analysis would be rather challenging because it requires 
one to solve the Laplace equation within a domain with a 
complicated shape in order to obtain the equilibrium con- 
figuration of the liquid-vapor interface. In the following 
we adopt this translational invariance and all extensive 
quantities are taken to be divided by the length of ex- 
tension in this direction. For example the expression for 
the energy of the system (see, e.f., Eq. ([9|)) denotes the 
energy per length. 

For our simplified model we introduce the additional 
approximation that the top of the plates is always at a 



vertical distance H from the substrate. This approxi- 
mation is justified if the height H of the plates is much 
larger than the distance L between vertical plates and if 
the clusters do not involve too many plates (see Fig.[3|). 

In order to calculate the adsorption we resort to a 
macroscopic model, in which the liquid-vapor interface 
is a section of a cylinder with its radius of curvature R 
given by Laplace's equation: 



R 



ApA/i 



(2) 



(A more sophisticated model would include a thin layer 
of liquid adsorbed on the walls [33 | . In a first approxima- 
tion we neglect this.) Our macroscopic (thermodynamic) 
approach implies that, for the non-tilted plates, only at 
capillary condensation the liquid meniscus forms a part 
of a circle, the radius of which is given by Eq. ^ and 
which meets the vertical wall with the contact angle 9. 
Actually, this is precisely what signals the occurrence of 
capillary condensation. 

The free energy of the system for a given configuration 
consists of a volume term (capturing that the system is off 
liquid- vapor coexistence), surface terms, and the elastic 
energy of tilting the plates. In principle this free energy 
has to be minimized in order to obtain the equilibrium 
configuration of the system and, from this, the adsorp- 
tion isotherms. However, the implementation of this pro- 
cedure is rather complicated. Moreover, for these types 
of systems it is not obvious that the y a re always able 
to reach the free energy minimum [4^ |45| | . Therefore we 
use instead a more ad-hoc approach, in analogy with Ref. 

Our aim is to determine the values of the chemical po- 
tential for which the clustering of the plates occurs. In 
between these values for the clustering events we treat 
the substrate as to be rigid and describe the adsorption 
without considering elasticity. To be specific, we cal- 
culate the values Aii^^\ Afi'^^\ AiJ,^^\ ... of the chemical 
potential for which the first, second, third, and so on clus- 
tering event occurs. The adsorption between these values 
of the chemical potential is taken to be the same as the 
adsorption on a rigid substrate with the plates fixed to be 
clustered correspondingly, i.e., two-by-two, four-by-four, 
etc. 

The first step is to determine if the plates cluster at all. 
According to the Kelvin equation 0] capillary condensa- 
tion between two parallel plates occurs at Afj, = Afi'^'^ 
(see Eq. [T|) . This equation is valid for condensation be- 
tween two infinite parallel plates. For a finite pore, i.e., 
for a capped capillary, the condensation transition still 
occurs for the same undersaturation A/i'^'^ but the tran- 
sition changes from first-order to continuous (albeit with 
a very steep increase of the adsorption) [s^l . Once again, 
the accurate description of this transition requires a more 
detailed model but this is of no relevance for the present 
analysis. In the following we consider thin plates, the 
thickness of which is much smaller than L which allows 
us to neglect it. However, capillary condensation does 
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not occur in all pores simultaneously because as soon as 
one pore fills, the surface tension pulls its confining plates 
together, causing the neighboring pores to open and thus 
preventing condensation in these. Therefore, in order to 
determine if the plates cluster, one has to balance the 
elastic, the volume, and the surface tension energies be- 
tween two individual plates filled by the liquid up to the 
top: 



E ^ 2fc(7r/2 



ApA^iH{L-x)+aiL-2x). (3) 



Here k is the elastic constant, a is the angle the plate 
forms with the horizontal, and x is the deviation of the 
top of the plate from the reference configuration a = 7r/2; 
a; > (x < 0) refers to closing (opening) the gap. Equa- 
tion dS]) assumes a planar instead of a bent liquid-vapor 
interface at the top of the pore, which represents a mini- 
mal requirement for the collapse of the plates, i.e., if the 
pores are not completely filled up to the rim the area of 
the liquid-vapor interface is larger. Therefore the config- 
uration of a flat interface at the top is the one with the 
largest possible, bending induced, free energy reduction 
associated with the liquid-vapor interface. This implies 
that if the surface tension cannot cause the bending of 
the plates even when the interface reaches up to the top 
everywhere, it never will lead to bending. 

Defining Sa = 7r/2 — a and using Sa ~ tan((5a) ^ x/H 
and solving for dE/dx = we obtain x* = + 

^^ff^. If X* < L/2 the plates bend due to the capiUary 
forces but do not touch each other and, as the chemical 
potential is increased further (i.e., is decreased), the 
neighboring pores with a wide opening fill without any 
further clustering. On the other hand, if the formal solu- 
tion X* of the equation dE/dx = fulfills x* > L/2, it is 
outside of the physically allowed domain for x* ; in fact 
this means that the plates would have to "cross" each 
other. Thus the condition x* > L/2 tells that the energy 
is minimized if the plates bend until they touch each 
other, at which point a hard-body repulsion prevents the 
"crossing" of plates. The above dependence of x* on A/i 
deserves a comment, as it might appear to be surprising 
that a larger A/j, leads to a larger x* . One must keep in 
mind that this relationship is valid under the assumption 
that the liquid-vapor interface is stuck at the top of the 
plates. The increase of A/x, while forcing the interface to 
stay at the top of the plates, causes the plates to bend as 
much as possible in order to reduce the volume of liquid 
as it is required for large undersaturations. However, this 
is not important for our analysis which is focused on the 
behavior of x* only near capillary condensation. There- 
fore the necessary condition for two individual plates to 
cluster is 



L < 



2k 



(4) 



After this first clustering, the system is formed by 
capped wedges. If we assume that in this first step (i = 1) 



the plates cluster pairwise, then by four {i = 2), followed 
by eight {i = 3) plates, and so on, the outermost tilt 
angle after the clustering of the i'^ order is 



tanai 



2H 



(2* - 1)L' 



(5) 



Note that a^^o = f ■ 

In the following we first consider a system of fixed clus- 
ters of order i. In order to obtain the corresponding ex- 
cess adsorption Ti we simply add the volumes Fc*'' and 
r;;' of liquid in a cluster of order i and in the remaining 

(i) 

wedges, respectively, multiplied by the total number TVc 
of clusters of order i in the system with TV plates where 



(0 



N/2^ after the i collapse. We have 



tan a," 



(6) 



which, within the present simple model, is independent 
of A/i. Thus we consider the situation that fluid trapped 
within a pore does not evaporate if A/i is increased. (To 
a certain extent this lack of evaporation upon increas- 
ing Afi mimics hysteresis.) A more refined model would 
take this evaporation into account, but this has no conse- 
quence for our main results. For Fw'' we have (see Eq. ([2]) 
and Eq. ([Tf)) in Appendix E| 



(ApA/i)2 



sin(ai 



I cosi 



cos ai 



The expression for Fw"* vanishes for A/i Ap,\ 

This signals 



(a, - 9) 



£2 

- — tan tti 



(7) 



4a^ 

(ApL)^tanQi 



sin(cti —8) cos 9 



(a. - 0) 



that for A/i > A/if^-* there is only a microscopically 
thin adsorbed film which is not captured by the present 
macroscopic description. Accordingly the expression in 



Eq. dZl) is augmented by r\i' (Ap > Ap]l') = 0. 

A more convenient quantity is the adsorption per unit 
area, i.e., per length if one divides by the extension of 
the system in the translationally invariant direction: 



— = (fW +F«^ — 
NL ^ ' NL 



T. 



2'L 



(8) 



Now we consider flexible plates. The clustering events 
occur if a configuration becomes linearly unstable. In 
order to calculate the value A/i*^*) of the chemical poten- 
tial at which a given configuration (i) leads to cluster- 
ing of order {i + 1) we check the stability of the system 
against a small perturbation which tilts the clusters an 
angle 6a towards each other. We consider the situation 
that the clusters tilt in pairs, so that the total energy 

i5tot((5a) = [Ei{ai + Sa) + Ei{ai - 5a)\ is given by 
(see Eq. (jAlli in Appendix \^ 
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ApAfj. 



C0S9 r /,x COS 

■ siii(Q;i + oa — u) — 



cos(ai + Sa) 



tanftti + Sa) 

4 a 



(ApA^)2 L 



cos(aj; — Sa) 
(ApAAi)2 



sin(ai — Sa — 9) 



— — tanfofi — Sa) 
4 cr 



(9) 



cos 9 ApA^ cos 9 ApAp 

7 ^ — ^ ~\~ J^' 



cos(ai + Sa) a 



cos{ai — Sa) a 



+ 2'+iA:(fo)2 + C«(a,), 



where the last contribution contains terms which do not 
depend on Sa and thus can be dropped upon varying Sa. 

For each i one has to check if the configuration Sa — Q 
is hnearly stable. By symmetry = is always an 
extremum of E\^x_. If it becomes linearly unstable the 
plates form clusters of order i + 1. This loss of stability 
can occur smoothly, with the plates tilting gradually to- 
wards each other, or abruptly, with the plates clustering 
instantaneously. Numerically we have always observed 
an abrupt clustering of the plates. 

It is too difficult to analytically determine the point at 
which the configuration with = is no longer linearly 



stable, i.e. 



dSa-' 



0. However, for given values of 



Sa=0 



the parameters we can proceed numerically to determine 
both the values A^^'^ of the chemical potential for the 
clustering events and the isotherms. 

Clustering and adsorption can proceed as long as the 
wedges between clusters are not filled by liquid. For a 
fixed cluster, this ceases ii h — h' = H (sec, c. f.. Fig. [71), 
i.e., at 



Ap = Ap* = 



a sin(Q;i — 



Ap(i7/tanQ;, + L/2) 



(10) 



with ai given by Eq. ([5]). For A^ < Ap* the liquid com- 
pletely covers the plates. In this case the surface tension 
does no longer promote the clustering of the plates be- 
cause the liquid-vapor interface is above the plates. Ac- 
cordingly the elastic forces will break the clusters apart, 
unless other forces, such as van der Waals forces, keep 
them together. 

In is also worthwhile to point out that the above de- 
scription automatically takes into account the possibility 
that, beyond a certain number of clustering events, the 
clusters are too rigid to be bent by the surface tension. 
In such a case the geometric configuration never becomes 
unstable and the liquid adsorbs until it completely covers 
the plates, as discussed at the end of Section HIl 



IV. DISCUSSION 

In order to translate the previous calculations into spe- 
cific results for the adsorption isotherms it is convenient 
to introduce dimensionless quantities, which we denote 



by an overbar such as E. To this end we measure the 
lengths in units of L, the energy per length in units of 
aL (the translational invariance of the system has been 
already taken into account) and the chemical potential 



in units of 



ApL ■ 



The elastic constant is also measured 



in units of aL and the adsorption is measured in units 
of L^. With this choice of units we are left with the pa- 
rameters H = H/L, which defines the geometry, 9, and 
k = which is a ratio of elastic and surface energies. 
It is not easy to estimate k for actual materials as we 
are assuming the plates to tilt around their base. In any 
case a large value of k signals a rigid substrate, for which 
the surface tension is not capable of bending the plates. 
On the other extreme, a soft substrate is described by 
a small value of k such that the surface tension forces 
cause extensive clustering. The spectrum of actual ma- 
terials spans the entire range of values, as can be inferred 
from experimental results, according to which clustering 
can either occur or not [i-Illlll,!!!]. It is to be expected 
that the interplay of elasticity and surface tension will be 
most pronounced for intermediate values of k. Accord- 
ingly, in Fig. [5] we present our results for the adsorption 
isotherms (in dimensionless units as described above) for 
k — 50, 9 = 0, and H = 100. For these parameters the 
necessary condition for clustering of two individual plates 

to occur, 1 < — — h Ap-^ (see Eq. HJ, is fulfilled for all 



2k 

(100)' _L A^ilM! 



50 



100 



k 

Ap, because 1 < 

Fig. [5] shows the adsorption isotherms for the first few 
fixed cluster configurations {i = 1,...,6) with Ti/N = 



2-^0'' +r{^'') where F^'^ and V'^' are given by Eqs. ^- 
([7]). In line with the discussion after Eq. ([7]) these 

isotherms are flat for Ap > Ap"^ with plateau values 

^^tttH — !• H/2 as i — > oo. For decreasing undersat- 
uration the isotherms increase up to Ap*{9 = 0) = 

2^-'e.(l+ef )^/^ 6 - 2?^, (see Eqs. © and (HU])) where 
the wedges become completely filled up and the cluster 
configurations dissolve; Ap*^^{9 = 0) = A~^H — > 0. In 
Fig. [5] this latter feature occurs outside the field of view 
(e.g., Apl « 1 and Vi{Ap{)/N > 60). When the clusters 
dissolve the adsorption jumps to the same value as for the 
rigid substrate, as sketched in Fig. 21 because elasticity 
no longer plays a role. 

If the cluster configurations are not fixed but flexi- 
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blc, at certain undcrsaturations A/i^*^ (determined by 

= 0, see 



d5o 



the linear stability analysis criterion 
Eq. ([T0|) above) the substrate configuration undergoes 
a first-order transition between clusters of order i (i.e., 
containing 2* plates) and clusters of order i + 1 upon de- 
creasing A/2, with corresponding jumps in the adsorption 
(see the black line in Fig. [S]). In between two jumps the 
adsorption isotherm of the fiexible substrate follows the 
isotherm of the corresponding fixed cluster. This leads 
to the step-like character of the black line in Fig. [S] 

Due to the choice 9 = capillary condensation occurs 
at AfL'^'^ = 2 (Eq. (jT])). Thus for the present system clus- 
tering sets in for this value of the chemical potential be- 
cause, as stated above, Eq. (|4]) is fulfilled. Therefore one 
has A/i'^°' = A/i^"^. The significance of elasticity is high- 
lighted by a comparison of the black line in Fig.[5]with the 
adsorption isotherm expected for the corresponding rigid 
substrate featuring the same plates but fixed in upright 
position as sketched in Fig. 2] At capillary condensation 
(A/2 = Afl'^'^ = 2) the adsorption for the rigid substrate 
jumps suddenly from a small value to T/N = 100 which 
is outside the field of view of Fig. [51 corresponding to 
the complete filling of the space between the plates. Our 
macroscopic model does not take into account the micro- 
scopic details of the fiuid-solid and fiuid-fluid interactions 
and therefore it cannot describe the precise shape of this 
isotherm, only the position of the jump at capillary con- 
densation, predicting for the present choice of parameters 
constant values of and 100 for the adsorption below and 
above the capillary condensation. A more detailed model 
(see, e.g., Refs. [3Jt33) is required in order to be able 
to determine the full shape of the isotherm. 

The isotherms for elastic substrates differ qualitatively 
from the isotherms of substrates composed of rigid clus- 
ters of fixed order i with a tilt angle given by Eq. (O 
(see also Fig. S]). As expected, the adsorption isotherm 
of the elastic surface exhibits abrupt jumps at each clus- 
tering event (i.e., at A/i = A/^*^*)), reflecting the change 
of the substrate geometry. The first two clustering events 
occur at undersaturations A/i very close to each other, 
because the large ratio H/L = 100 implies that the an- 
gles Qfi of the wedges formed in the first clustering events 
are very close to 90° and thus in Fig. [S] the differences 
between these values of A/x*^*^ are barely visible. Another 
notable feature is that for decreasing A/i the adsorp- 
tion is not necessarily monotonically increasing. This 
counter-intuitive behavior is due to liquid evaporating 
(as indicated by the downward arrows in Fig. ^ from 
the capped wedges, which form when a cluster collapses 
upon decreasing A/i. An experiment in which the chem- 
ical potential is increased slowly, in order to facilitate 
thermal equilibrium and evaporation of liquid into the 
vapor phase, is expected to pick up this peculiarity. 

It is also informative to monitor the filling height h — 
h — R{1 — cos{a — 9)) of the liquid in the center of the 
capped wedges (see Fig. [7] in Appendix A) as function 
of the order i of clustering (i.e., at A//^*)). For i = 1 to 
5 the heights in units of L are 0.001,0.73,2.31,7.1, and 



34.2 compared with the maximal value H/L = 100. The 
first few values are rather small, reflecting the fact that 

= 50 describes rather flexible substrates which cluster 
easily. For these low orders, the plates cluster almost as 
soon as an interface is formed. 

As a final comment we point out that our choice of 
fc = 50 was not accidental. For a smaller value of k the 
substrate is so flexible that clustering occurs even before 
the interface spans the space between two plates. The 
formation of the interface associated with the liquid ad- 
sorbed in the corner formed at the point of contact of 
the plates with the substrate is sufficient to cause clus- 
tering. A signature of this feature is that the filling height 
h~h — R{1 — cos(a — 9)) discussed in the previous para- 
graph becomes negative if one applies the approach de- 
scribed before. This shortcoming for small k could be 
fixed easily by taking into account the adsorption in the 
corners formed between the plates and the underlying 
substrate (see Fig. [5]). However, this extra step would 
not add any qualitatively new feature and would compli- 
cate our simple analysis unnecessarily. 



V. CONCLUSIONS 

Our analysis shows how a responsive substrate and the 
coupling between wetting and substrate geometry can 
lead to new phenomena. For the specific system consid- 
ered here the adsorption isotherms exhibit steps at un- 
dersaturations for which the substrate changes geometry 
due to the adsorption of liquid. 

For these phenomena adsorption isotherms are partic- 
ularly revealing in that they are directly accessible to ex- 
periments, as already proven for geometrically patterned 
substrates |47| . and they provide an expressive charac- 
terization of the wetting properties of substrates. Our 
choice of the substrate might not be the easiest to be 
prepared experimentally but in view of grea t advances in 
the fabrication of patterned substrates [43, it might be 
possible to fabricate such a substrate or a similar one. A 
substrate with fioppy plates could be created more easily, 
but the corresponding theoretical analysis will be more 
difficult, because in this case the bent plates will not form 
straight wedges but rather more complicated shapes. 

We have used a simple macroscopic model in order to 
calculate the adsorption. It is likely that for a direct 
quantitative comparison with possible experimental ad- 
sorption isotherms our simple model is insufficient. A 
more accurate model should incorporate more realistic 
details such as van der Waals interactions [l^] and more 
realistic elastic features of the patterned substrates. Such 
refined descriptions will be most successful if they are tai- 
lored to specific substrates used in actual experiments. 

We end with mentioning a prospective application of 
our results. So far we have considered "passive" sub- 
strates. The geometry changes with adsorption and in 
turn this change is determined by the amount of liquid 
adsorbed and by the characteristics of the substrate fixed 
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FIG. 5: (Color online) Adsorption isotherm of an elastic, patterned substrate (black line) and isotherms for substrates composed 
of rigid, clustered plates (dashed lines) with clusters of the order i = 1, ...,5 from right to left at the top, and from bottom 
up at the right. The first three configurations of the substrate are sketched, along with pointers indicating the Ap, intervals 
within which they occur. The clusters of order i containing 2* plates are stable for A/i'*' < Ap, < A/i''"^^ (with Ap'^^' = oo 
for i = 0). The coUapse of the clusters of order i takes place at Afi^°^ = 2, Ap'^^ = 1.990, A/i'^' = 1.928, A^i'^^ = 1.663, 
Ap,^'^^ = 0.901, and A^i'^^ = 0.148. The first two collapses occur very close to each other and consequently they are barely 
discernible in the isotherm. The systems are characterized by the dimensionless parameters L = 1, H = 100, k = 50, and 
9 = 0. See the main text for a detailed discussion. 




FIG. 6: (Color online) Liquid-vapor interface which does not 
span the space between two clusters but which connects both 
plates with the substrate surface. Our model does not take 
into account the adsorption in the corners, but this is impor- 
tant for extremely fiexible substrates. 



cities. Such controllable surfaces might find applications 
in, e.g., "lab-on-a-chip" devices. 
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via its fabrication. One can instead imagine a more "ac- 
tive" substrate by using, e.g., electric fields in order to 
control how the substrate structures collapse. Indeed 
electric [i^ and magnetic fields jHO] have already been 
used for control of the collapse of elastic substrates. In 
this setup one could change externally the effective stiff- 
ness of the substrate, and thus control its wetting prop- 



Appendix A: Expression for the free energy 

Within our model the calculation of the free energy for 
a given geometric configuration and chemical potential 
reduces to geometrical considerations. The free energy is 
the sum of volume, surface, and elastic contributions: 



E = ApA/i X (volume of liquid) + cr x (liquid- vapor surface area) 
-f- ais X (liquid-solid surface area) + a^s x (vapor-solid surface area) 

-I- elastic energy. (Al) 
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FIG. 7: Geometry of the problem depicting the relevant 
lengths and angles. The gray area denotes the liquid. The 
actual surfaces are indicated by thick black lines. One has 
I = i?sin(a - 0), h! = M|£L^^ and h = Ztana. Note that 
for a given geometry (i.e., q and L) and given materials (Q) 
the filling height h — h' ~ R{1 — cos(a — 6)) of the liquid in 
the center-line of the pore is determined by the radius R of 
curvature, which in turn is given by the undersaturation Afj, 
(see Eq. Q). 



Equation (|A1I) assumes that the bottom of the cavity 
and the plates consist of the same solid material so that 
both give rise to the same surface tension <t\s. Fig. [7] 
illustrates the geometry of the problem and defines all 
relevant lengths and angles. In the following we omit 
the linear extension of the system in the translationally 
invariant direction. 

Within a harmonic approximation the elastic energy of 
a single bent plate forming an angle a with the horizontal 
substrate is 

E,i = k{7r/2-af, (A2) 
so that the elastic energy of a cluster of the order i is 

i;«=2fc^(V2~a,f. (A3) 



Accordingly, the additional elastic energy required for 
bending a cluster of the order i whole is 



2-1 



-E 



= 2fc^(<5a, 



(A4) 



The bending of a cluster as a whole leads to a displace- 
ment of the top of the bundle, characterized by the cor- 
responding angular deviations 5aj of the participating 
plates. As an approximation we take 5aj = 5a, which is 
valid if the lateral displacement of the top of the plates 
due to bending is small compared with the height H of 
the plates. This leads to 



5E^^ = 2'k{da) 



(A5) 



Concerning the contribution to the volume terms of 
the free energy one infers from Fig. [8] that the volume of 
liquid (divided by the linear extension of the system in 
the translationally invariant direction) in a wedge after 
the i^^ clustering is 



r(i) 



lh-h'L/2 

- [{ai - e)R^ - i?2 sin(a 
sin(Q;i — 9) cos 



R 



cos ai 



6) cos(q; 
{a^ - 9) 



0)lM) 



L2 

tan ai 

4 



Likewise the surface contributions follow from inspection 
of Fig. [71 The surface area of the liquid-vapor interface 



FIG. 8: Representation of the volume of liquid adsorbed. 



is 

S'^ = 2{a, - e)R. 
The surface area of the liquid-solid interface is 
2 



cos a. 



■{Rsm{ai-e) -L/2). 



(A7) 



(A8) 



Finally, the surface area of the vapor-solid interface is 
2H 



sm ai 



L. 



(A9) 
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Adding these terms and using Laplace's equation one obtains for the total free energy per number of 
(Eq. ^) and Young's equation clusters 

cosfl^^^^l^ (AlO) 



cos a, 



■ sin(Q!i 



L2 (ApA/i)2 



tan a,- 



cos 9 ApA/i 



cos a, 



(All) 



In Eq. lAllI we have omitted contributions which are in- 
dependent of a," and the contribution ?^ to the vapor- 

^ sin ai ^ 

solid surface energy. This latter term is an artifact gener- 
ated by assuming a constant value of H for the position 
of the top of the plates. Moreover we have dropped the 
contribution vi^^ ApAfx from the above expression. These 
terms might seem to be relevant as they depend on ai. 



However, the important dependence resides in l/cosai 
because we always consider the values of ai to be close 
to 7r/2 and thus the terms we have discarded are quan- 
titatively unimportant. In the same spirit one could also 
drop the term proportional to ai — 9, which will be can- 
celed in further calculations anyway. 
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